Fermi-pocket-vanishing phase transition in the ferromagnetic phase of the Kondo lattice model 
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By performing extensive dynamical mean-field theory calculations with the numerical renormalization group 
as the impurity solver, we establish the low-temperature phase diagrams of the spin-1/2 and spin-1 Kondo 
lattice models as a function of the conduction-band filling and the exchange coupling strength in the regime of 
ferromagnetic RKKY interactions (i.e., for low and moderate electron density). We show that both models have 
two distinct ferromagnetic phases separated by a continuous Lifshitz transition of the Fermi pocket vanishing 
type: one phase has a true gap in the minority band, the other only a pseudogap. The two phases can be 
experimentally distinguished by their magnetization curves; only the first phase exhibits magnetization rigidity. 
We find that, quite generically, ferromagnetism and Kondo screening coexist rather than compete, both in spin- 
1/2 and spin-1 models. We compute the Curie temperatures and establish a "ferromagnetic Doniach diagram" 
for both models. 
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Materials with competing interactions have complex low- 
temperature phase diagrams with quantum phase transitions 
separating the different ground states Heavy-fermion 
compounds containing ions with partially filled / shells show 
a particularly wide range of ordered phases llsj-^l. The 
paradigmatic model for these materials is the Kondo lattice 
model (KLM) OhSO which describes a conduction band of 
itinerant non-interacting electrons and a lattice of local mo- 
ments; the itinerant and localized electrons are coupled at each 
site by an antiferromagnetic exchange interaction J. For large 
J, the itinerant electrons screen the local moments in a lattice 
version of the Kondo effect. The resulting state is paramag- 
netic with Fermi liquid properties, albeit with strongly renor- 
malized parameters. For small J, however, the conduction- 
band electrons behave as carriers of the long-range RKKY 
interaction which leads to magnetic ordering of the local mo- 
ments. The two regimes are separated by a quantum phase 
transition at critical ,/*. This competition between the lo- 
cal Kondo effect and the intersite exchange interaction is de- 
scribed by the Doniach diagram ifioll . While the Kondo tem- 
perature is an exponentially increasing function of J, Tk oc 
exp(— 1/ pj), the Neel temperature increases at first quadrat- 
ically with J, but then it peaks and decreases to zero at J* as 
the strong Kondo screening takes over. The simplest version 
of the KLM with spin-1/2 local moments indeed has an an- 
tiferromagnetic (AFM) ground state (Neel order) for small J 
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m, UCo^^Sbz im, NpNiSiz 123|], 
In addition, there are 



near half-filling of the conduction band llllLll2ll . The nature 



of the quantum phase transition at the critical J* has been in- 
vestigated using a variety of methods, the most accurate of 
which confirm that the transition is second order (quantum 
critical) and indicate that it involves a change of the Fermi sur- 
face topology II13I - I15II . In the spin-1 KLM, there is no phase 
transition at half-filling and the AFM phase extends to large 
values of J. 

While cerium compounds commonly have antiferromag- 
netic ground state, some are ferromagnetic (FM): CeRu2Ge2 
ifiell . Celn2 USE!], and CeRu2Al2B iH]. A number of ura- 
nium and neptunium heavy-fermion materials are also FM: 



strong indications of robust coexistence of the Kondo ef- 
fect and ferromagnetism, in particular in U compounds. In 
Refs. JzsUigll it has been proposed that an appropriate mini- 
mal model for this behavior is the spin-1 version of the KLM, 
where the conduction-band electrons underscreen the local 
moments, while the residual moments order ferromagneti- 
cally. FM order appears for low and moderate electron fill- 
ing in the conduction band Mean-field analysis 
predicts two phases: for small J the stable phase is a FM reg- 
ular metal (i.e., the hybridization mean-fields are zero, while 
the magnetization mean-fields are finite), while for large J 
there is a transition to a FM heavy metal (i.e., both mean-fields 
are non-zero). Dynamical mean-field theory (DMFT) calcula- 
tions demonstrated that the standard spin-1/2 KLM also has a 
FM order coexisting with (incomplete) Kondo screening [34 1 . 
Furthermore, this phase turned out to be a half-metal (the mi- 
nority spin band is gapped Is?]) and there is a commensura- 
bility condition relating the magnetization to the band filling 
1 3411 . which can be interpreted within the hybridization mean- 
field picture as being due to completely filled minority-spin 
lower band |3^ 37 1- A recent mean-field analysis of the spin- 
1/2 model suggested the presence of several different ferro- 
magnetic phases 1380. So far, however, a single FM phase has 



been identified in the DMFT calculations ll32Ll33[l 



These findings open a number of important questions: 
What is the relationship between ferromagnetism and Kondo 
screening: do they compete or coexist? If there is some degree 
of competition, how does it manifest? What is the minimal 
model for studying these effects, spin-1/2 or spin-1 KLM? Is 
there a quantum phase transition between different FM states 
also in the spin-1/2 model? What is the nature of these tran- 
sitions and what are their experimental signatures? And, fi- 
nally, which aspects of the static mean-field analysis are cor- 
rect and which must be revised in more accurate dynamical 
treatment? To answer these questions we have performed ex- 
tensive DMFT 1I39I] calculations using the numerical renor- 
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Figure 1: (Color online) Phase diagrams of spin- 1/2 and spin-1 
Kondo lattice models for small and intermediate occupancies n. 
Phase A is a ferromagnetic half-metal phase with strong Kondo effect 
where the minority band is gapped. Phase B is an itinerant ferromag- 
netic phase with a pseudogap. Phase C for spin- 1/2 model indicates 
the region where the true ground state is expected to be a charge 
density wave. For veiy small occupancies, the calculations fail to 
converge due to limitations of the impurity solver. 



maUzation group (NRG) as the impurity solver 11401-145 1 
We consider the Kondo lattice model 
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which describes a single-orbital conduction band with disper- 
sion Lj = efc, and a lattice of local moments described by the 
spin-S* operators S^; is the conduction-band spin-density at 
site i, and J is the antiferromagnetic Kondo exchange cou- 
pling (J > 0). We focus on the Bethe lattice that has a semi- 
circular density of states with bandwidth 2D. 

In Fig. [T] we present the main result of this work: the phase 
diagrams of the spin- 1/2 and spin-1 KLM as a function of 
the band filling, n < 0.5, and the strength of the exchange 
coupling, J. For both spins, 1/2 and 1 alike, we find two 
different ferromagnetic phases. The one at large J (phase A) 
corresponds to the ferromagnetic half-metal phase described 
by Peters et al. ll34ll . The corresponding spin-resolved spec- 
tral functions for the 5 = 1 model are shown in Fig.|2] panel 
A. The minority spin band is gapped ll34ll . while the majority 
band only exhibits a weak hybridization pseudo-gap charac- 
teristic of the Kondo lattice systems |4^ 47]. The phase at 
small J (phase B) is not gapped, but there is a pronounced 
pseudogap just below the Fermi level in the minority band, 
Fig El panel B. The spectral functions for the 5 = 1/2 model 
are qualitatively the same. The spectra thus suggest the occur- 
rence of a Lifshitz transition: there is no change in the sym- 
metry, but the Fermi surface of the minority band shrinks to 
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Figure 2: (Color online) Spin-resolved conduction-band local spec- 
tral functions Ac.cr for the spin-1 KLM in the ferromagnetic half- 
metal phase (A) and in the itinerant ferromagnetic phase (B). The ar- 
row indicates the main effect of decreasing interaction J: the lower 
edge of the upper hybridized band shifts to lower frequencies. The 
left insets in both panels show the /-level spectral functions Af,^ 
defined through the imaginary part of the scattering T matrix. The 
right inset in the upper panel shows the spectral functions in the full 
frequency interval. 



a point and disappears as one goes from phase B to A. We 
emphasize that the two phases exist both for spin- 1/2 and for 
spin-1 models and have similar properties; clearly, the value 
of the spin does not play a crucial role in the FM ordering. 

The transition point J* is a non-monotonic function of fill- 
ing that peaks at n ~ 0.2 and n ^ 0.25, respectively. Near 
n ^ 0.4 we observe a change of behavior in the small- J 
phase. For 5=1/2 KLM, this is the parameter regime where 
a charge-density wave (CDW) ground state is predicted 113 311 . 
but not allowed for in our calculations. 

In order to study the nature of the phase transition between 
A and B more carefully, we plot in Fig. |3]the evolution of 
the magnetization (both total, and for the conduction-band c 
and localized / electrons separately) and of the quasiparticle 
renormalization factors 



Z„ = 



(2) 



as a function J across the transition point J*. The frozen 
magnetization in phase A is given by a generalization of the 




Figure 3: (Color online) Total, conduction-band c-level and localized 
/-level magnetizations (top panels) and the spin-dependent quasi- 
particle renormalization factors Za (bottom panels) across the phase 
transition, indicated by the vertical dashed lines. The magnetiza- 
tion is here defined as the expectation value of the spin operator 
without the —guB factor: mj = (Sz), rric — (n-|- — ni)/2, 
rritotai = m/ + rric. In the plots, mtotai and m/ are shifted by 
ms defined in Eq. ([3}. 



spiii-1/2 KLM result from Refs. SIMS: 



ms = {2S -n)/2. 



(3) 



The magnetization is continuous, but there is a change of slope 
in TO/, while rric has a continuous first derivative. This is in 
disagreement with the static mean-field analysis for S = 1, 
which predicts a large jump in the magnetization ll27ll . The 
renormalization factors Za- for both spin orientations are con- 
tinuous and finite across the transition (note that in the mi- 
nority band of phase A there are no quasiparticles, but 
can formally still be defined). There is no criticality in this 
spin selective metal-insulator transition. From these results 
we may conclude that this is a continuous Lifshitz transition of 
the Fermi pocket vanishing type SmiiMIll. The change 
of the Fermi surface topology is continuous in the sense that 
there is no sudden reorganization of the Fermi surface. Deep 
in the phase A, the majority electrons become weakly corre- 
lated {Z saturates at a value of order 0.5). 

For very large J, there is another Lifshitz transition to a 
non-gapped phase 115 211 that we denote as B'. While in the BA 
transition, the chemical potential is located at the bottom of 
the upper hybridized band, in the AB' transition the chemi- 
cal potential is located at the top of the lower hybridized band 
at the transition point. In other words, while BA corresponds 
to the vanishing of electron pocket, AB' corresponds to the 
emergence of hole pocket. The AB' transition occurs at un- 
physically large values of J in excess of the bandwidth (the 
Schrieffer- Wolff transformation mapping the periodic Ander- 



Figure 4: (Color online) Magnetization in longitudinal external mag- 
netic field. The dashed lines indicate the value of the frozen mag- 
netization ms- The ,g-factors are assumed equal for c and / levels. 



son model to the KLM ceases to be valid), therefore we do 
not plot them in the phase diagram. For spin- 1/2 model with 
n = 0.3, this AB' transition occurs at J/D = 2.25. For even 
larger J, the system eventually becomes paramagnetic in the 
ground state. For n = 0.3, this occurs for J/D ^ 3.4. 

The static mean-field also predicts four distinct phases at 
T = t38il : ferromagnetic, "spin-polarized" Kondo, "spin- 
non-polarized" Kondo, and paramagnetic phase, the first three 
roughly corresponding to our B, A, and B' . We find, however, 
that all four phases that we find exhibit the Kondo effect, i.e., 
our phase B is not pure ferromagnetic. Furthermore, the Lif- 
shitz transitions we find are all continuous: there are no jumps 
in any of our results. Finally, we note that the DMFT predicts 
that deep inside B and B' phases there are pseudo-gaps rather 
than gaps (this is due to non-zero imaginary part of the self- 
energy, i.e., due to correlation effects). The most surprising 
outcome of the DMFT calculations is, therefore, the gradual 
emergence of true gaps from pseudo-gaps as the gapped phase 
A is approached from B or from B', while the mean-field re- 
sults are closer to the rigid-band picture. 

Does the existence of multiple phases indicate a compe- 
tition between the RKKY interaction and the Kondo effect? 
Some degree of antagonism is suggested by the fact that the 
/-shell magnetization m / has a minimum at the BA Lifshitz 
point where both tendencies are expected to be equally strong 
and, furthermore, it could be argued that m / increases with J 
in phase A only because Kondo screening is rendered incom- 
plete by the opening and widening of the gap. Nevertheless, 
this competition does not imply mutual exclusion and most 
results rather support the notion of robust coexistence. 

Experimentally the phases A and B can be distinguished 
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Figure 5; (Color online) Temperature dependence of the magneti- 
zation, resistivity and heat capacity for the spin- 1/2 Kondo lattice 
model in phases A and B. The horizontal axis is rescaled by the 
Curie temperature Tc- Resistivity is in units of po = 27re^<l>(0) /hD, 
where "1? is the transport integral. Heat capacity curve was obtained 
by differentiating a piecewise interpolation of the numerical results 
for the total energy. 



by their magnetization curves. In phase A, mtotai remains 
pinned to 7715 for a finite range of the field strength, while in 
phase B the susceptibility dM /dB near zero field is finite, see 
Fig. |4] For sufficiently strong field, a gap opens in the mi- 
nority band in phase B, too. This effect can be understood 
within a rigid-band picture which holds to a first approxima- 
tion. For very strong field, the magnetization is reoriented in 
a first-order spin-flop transition which preempts another Lif- 
shitz transition. 

In Fig. |5]we plot the temperature dependence of key ther- 
modynamic and transport properties in phases A and B. We 
find that the magnetization in phase B remains essentially 
pinned at ms until T becomes of the order of the gap, while 
it has a finite temperature-derivative at T = in phase A. 
This difference is, however, small. The resistance p increases 
in both phases up to the Curie temperature Tc, then it de- 
creases approximately as a power-law T~'^-^, not logarithmi- 
cally. The heat capacity c has a jump discontinuity at Tc- Sim- 
ilar features are indeed observed experimentally, for example 
in Refs. iH |22|,|23t], although the simple KLM does not cap- 
ture the full complexity of real materials. 

We summarize the behavior of both Kondo lattice models 
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Figure 6: (Color online) "Ferromagnetic Doniach diagram" for spin- 
1/2 and spin-1 Kondo lattice models. The arrows indicate the posi- 
tion of the Lifshitz transitions. 



in the form of a "ferromagnetic Doniach diagram" in Fig. |6] 
We plot the single-impurity Kondo temperature (which does 
not depend on the impurity spin 115 311 ) and the Curie temper- 
ature Tc for each model. The Curie temperature has no ob- 
servable feature at the Lifshitz transition points (indicated by 
the arrows) and in this parameter range it does not go to zero 
nor even decrease, unlike the Neel temperature in the anti- 
ferromagnetic case. These results indicate that the Kondo ef- 
fect indeed favors FM ordering through the RKKY interaction 
|26l|33l|34t]. Apart from the (approximately) factor of two dif- 
ference, there is no pronounced difference in Tq of spin- 1/2 
and spin-1 models. 

We conclude by answering the questions raised in the in- 
troduction. The ferromagnetism and local Kondo singlet for- 
mation do not compete in the same way as the antiferromag- 
netic ordering and the Kondo effect do; there is no Kondo 
breakdown and no criticality, but rather a continuous filling 
of the lower minority band and the disappearance of the elec- 
tron pockets. We thus find robust coexistence of FM order 
and Kondo screening in both phases. This is the case for both 
spin- 1/2 and spin-1 Kondo lattice models, thus the physics 
of Kondo underscreening does not need to be invoked to ex- 
plain the magnetic ordering. Both models have qualitatively 
the same phase diagram in this range of filling and exchange 
coupling strength. The Lifshitz transitions are observable in 
the temperature and magnetic-field dependence of the mag- 
netization. The static mean-field appears to be valid at the 
qualitative level, however to properly describe the real nature 
of ferromagnetic phases and transitions it is necessary to take 
into account dynamic effects, as in the DMFT treatment. 
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